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The rise and the fall of the QS population are explained in terms of massive black 
holes forming/accreting during the assemblage of the host galaxies, and of accretion 
rekindled by interactions of the host with companions in a group. We compute the 
LFs out to 2 = 6. We also predict the masses of relict black holes to be found in 
many galaxies. We compare the histories of the QSO and of the star light. 



1 The rise and fall of the quasars 

The body of radio, optical and X-ray evidence supporting a genuine, if broad 
peak in the number of bright QSs around z ~ 3 has been reviewed by Shaver 
ct al. 1996. The QS population is found to rise in number during the first 
few Gyrs of the Universe life, down to epochs z ~ 3; in the QS jargon, this 
is a negative DE toward increasing z. Later than z 3 the population falls 
in cosmic epoch on a scale of ~ 2 — 3 Gyr. This course has been interpreted 
as dominated by LE, but at least two additions are called for (see La Franca 
& Cristiani 1997): at the faint end the LFs grow higher with z, and require 
a considerable positive DE; the bright end of the LFs is flatter at low z, with 
more numerous bright objects than previously recognized. 

That QSs as individual sources are powered by gas accretion onto a massive 
black hole has been argued many times, and direct evidence is now mounting 
(see Rees 1997). But accretion limited only by radiation pressure exponentiates 
on the short time scale ijtE — 4 10 7 yr. The coherent rise and fall of the QS 
population over most of the Universe lifetime requires extended coordination. 
We show this to be provided by the hierachical growth of structures, first 
as the host galaxies are built up, and then as groups form where the hosts 
interact with companions. The transition takes place just at z ~ 3; but in 
both dynamical regimes the symmetry of the gravitational potential is broken, 
and this provides conditions for fueling or for refueling the BHs. 



2 BHs fueled by galaxy formation, and refueled by interactions 

The hierarchical cosmogony envisages the typical density perturbation collaps- 
ing and virializing at the epoch t to have a mass growing as M c oc t 4 /(™+ 3 ) 
when D, = 1; for Cold Dark Matter perturbations n ~ —2 applies. A similar 
trend is retained when Q a < 1 until the growth freezes out at 1 + z ~ l/f2 . 



1 



A subgalactic building block of DM with M ~ 10 9 M Q allows a 10 6 M BH 
to form, involving a baryon fraction e„ ~ 10~ 3 (HR 93). The galaxy assemblage 
goes on chaotically with angular momentum j 7^ const, the gas accretion is only 
Eddington-limited and produces luminosities L ~ Le oc Mbh ~ e*M. 

The transitional mass from a large galaxy to a small group is 5 10 12 M . 
In the hierarchical cosmogony this corresponds on average to 2, ~ 3 ± 0.5, 
depending on the cosmology and on the detailed perturbation model. Such 
values match the range where the bright QSs peak, corroborating our view. 

In groups, many simulations (see Governato et al. 1996) have shown in- 
teractions to be frequent and effective, due to the high density of galaxies 
n g cx pu{z), and to the velocity dispersion V cx M^ 1 ™)/ 12 ^ t being still 
close to the galaxian v g . Nearly grazing encounters occur on the time scale 
r r <~ l/n g nr g V, and produce merging (Cavaliere & Menci 1997) or other 
strong interactions which perturb the potential and again cause j 7^ const. In- 
deed, aimed simulations (see Barnes & Hernquist 1991) show the gas in both 
interaction partners to lose j and be driven down to the main galactic nucleus. 

This view is supported by the extensive astronomical evidence showing 
individual QSs in hosts either engaged in merging and in current interactions, 
or having close companions even submerged within the galaxian body (see 
Hutchings & Ncff 1992, Disney et al. 1995, Bahcall et al. 1996, Hasinger et al. 
1996). The statistics ranges from a lower limit ~ 15% (Rafanclli et al. 1995) 
up to ~ 1/2 of the hosts (Bahcall et al. 1996); the average QS environment 
comprises ^10 — 20 galaxies (Fisher et al. 1996). 

3 The evolution of the luminosity function 

To include such bimodal fueling of the same engines, our LFs (see also CPV 
97) comprise at any L, z two components: N(L, z) — N\(L, z) + N^L, z) . 

N\ accounts for BHs newly formed on the scale td yn of the host buildup a 
process dominant for M ^ 5 10 12 M Q and z ^ 3. It is computed on the basis 
of the Press & Schechter 1974 mass function N(M,z), using the Eddington- 
limited L cx M similarly to HR 93, but with two crucial differences: we let Mbh 
follow the hierarchical growth of the host galaxy; we consider the probability 
oc pff(z) of forming BHs as a prefactor of N(L, z) rather than of M BH oc L. 

N2 comprises BHs reactivated by interactions, and is dominant in struc- 
tures (especially in groups) with M > 5 10 12 M Q , at z 3. The BHs restart 
their bright career from a low luminosity distributed after f(L); the reactiva- 
tion of the N r dormant BHs occurs with probability oc l/r r . Now the activity 
is supply-limited, reaches a top Lb to statistically fade off on the time scale 



2 




6 7 8 91011 1 2 3 4 



log M/M s z 

Figure 1: The LFs at z = 6,5,4.5,3 (panel a, from bottom to top), and at z = 2.4,1,0.5 
(panel b, from top to bottom) computed for QS forming in the critical universe from tilted 
CDM density perturbations, and then accreting host gas upon interactions; optically selected 
data as in CPV 97. Panel c shows the mass distribution of dormant BHs expected from this 
model at the center of many normal galaxies; for the data, see CPV 97. Panel d compares 
the history of the light density from QSOs with bolomctric L > 10 45 erg/s, to that of the 
star light produced by Madau 1997; the normalization of the former is actually ~ 1/30. 

tl oc Lb/L. All that is expressed by the rate equation 

d t N 2 + d L (LN 2 ) = f(L) N r /r r - N 2 /t l . (1) 

To close the argument, iV r is made up by the last term integrated over L, 
and by the newly formed BHs described by Ni. Thus N 2 arises from Ni and 
requires no independent normalization; Ni is normalized to the data at z = 4, 
and implies today about one relict BH per few bright galaxies. 

4 Conclusions and discussion 

Fig. la shows the LFs at high z. The prediction at z = 6 is sensitive to the 
threshold S c for collapse in the Press & Schechter formula; this is taken here to 
decrease slowly with z from the value 1.7, as indicated by recent simulations. 
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Fig. lb represents the low-z behavior when the accreted gas is provided 
mainly by the host reservoir (say, a constant fraction used up in each inter- 
action); then -Mg as /M gas ~ t^ 1 ~ -L b /L b holds, to yield L b (t) oc t~ to l Tro . 
With r ro ~ 6 Gyr (scaled to groups from the classic census of local interact- 
ing galaxies by Toomre 1977), the result is L b oc (1 + z) 3 . This implies LE 
dominant at z < 3, and LFs flattened by L increasing over the flyby time 
2r g /V ~ few 10 8 yr. But fig. lb shows that quite some DE also occurs, since 
the frequency of the effective interactions peters off when groups aggregate 
into clusters with high V. Fig. lc shows the mass distribution predicted for 
the relict BHs in many local, currently inactive galactic nuclei. 

Alternatively, the gas may reach the nucleus of a large galaxy mainly from 
satellites cannibalized out of an initial retinue gradually depleted. Then L b ~ 
const obtains with no LE, while the DE is enhanced; the gas mass oc M sat 
follows the satellite distribution and yields steeper LFs at the faint end. 

In either case, one ends up with the QS light history (of gravitational ori- 
gin) shown in fig. Id, and compared with the star light history (of thermonu- 
clear origin) produced by Madau 1997; at z ~ 1 this is mainly contributed by 
the faint blue galaxies, which Cavaliere & Menci 1997 interpret as starbursts 
in dwarfs interacting in LSS. The similar run at such z corroborates the view 
that interactions in dense environments (condensing LSS or virialized groups, 
respectively) drive much of the light illuminating the z ~ 1 Universe. 
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